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ABSTRACT 


If a Pd alloy containing a substitutional base element (Al) is exposed to an 
oxidizing environment at high temperatures, internal oxidation of the base element may 
occur if the solute is present in small quantities. Internal oxidation of Pd-Al alloys results 
in alumina (AI2O3) precipitates within the Pd-Al matrix, this produces Pd-Al203 
composites. The shape of alumina precipitates varies as a function of temperature. The 
thermal expansion coefficient and Young’s Modulus of Pd matrix (a m = 11.5 x KT 6 K' 1 , 
Em = 123 GPa) is different from that of the AI2O3 precipitate (ctm ~ 8.0 x 10" 6 K" 1 , Em = 
380 GPa). Therefore residual stresses are generated in the composites upon cooling from 
the oxidation temperature due to difference of thermal expansion coefficients of Pd and 
AI2O3 precipitates. It is known that Hydrogen is absorbed by Pd-matrix and the presence 
of these residual tensile stresses in the Pd matrix results in solubility enhancements. This 
demonstrates the nature of stresses developed in the composite and therefore hydrogen 
has been used as a probe for estimating residual stresses in Pd-alumina composites. 

The strains and stresses developed in the Pd-Al203 composite have been 
thoroughly studied by Finite Element Method (FEM). In this study, the FEM has been 
applied to study the thermal strains and stresses at every point generated in the Pd-Al203 
composite. Two types of composites have been utilized to study the strains and stresses 
developed in the composite. The first model considers spherical AI2O3 precipitates 
distributed in an infinitely spherical matrix. This model has been utilized to obtain 
strains, stresses, thermal residual stress, solubility enhancement ratio and Helmholtz free 
energy. 

The second model considers cylindrical AI2O3 precipitates in an infinite 
cylindrical matrix. This model is a plain strain case. FEM has been applied for various 
shapes and positions of the precipitates. Strains, stresses, solubility enhancement ratio 
and Helmholtz free energy are obtained from this model. 

Strains obtained Pd matrix are always higher than AI2O3 precipitates due to high 
thermal expansion of Pd matrix. Higher thermal residual strains and stresses are 
developed in the composite with higher temperature difference. Stresses in AI2O3 



precipitates are higher than Pd matrix, this is due to higher young, s modulus of 
precipitates. Tangential stress for AI2O3 are tensile where as compressive stresses are 
developed in Pd matrix. Thermal residual stresses result in solubility ratio, higher ratios 
are obtained for higher temperature difference. Helmholtz free energies (HFEs) solely 
depend on strain and the temperature of the composite but not on the stresses. As the Pd 
strains are always higher than AI2O3 so their HFEs are always higher. 
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Chapter 1 


INTRODUCTION 

One method to obtain Pd-A^Oj composites is by internal oxidation of Pd-Al alloy 
at high temperature. Oxygen dissolves in the Pd-Al alloys and internally oxidation of A1 
produces alumina (AI 2 O 3 ) precipitate in the nearly pure Pd matrix. The thermal expansion 
coefficient of Pd matrix (a m = 1 1.5 x 10* 6 K' 1 , E m = 123 GPa) is different from that of the 
AI 2 O 3 precipitate (a m = 8.0 x IQ' 6 K' 1 , E m = 380 GPa). Thermal residual stresses 
developed when ideally oxidized alloys are cooled from high temperature to low 
temperature due to mismatch of material properties. These thermal residual stresses can 
be modeled theoretically also. 


1.1 Hydrogen Solubilities in Composites 


The stress distribution in the Pd-AbCb composite is complex. Hydrogen can be 
used as probe to understand the stress generated within the composite [1]. Hydrogen 
behavior in Pd is well-known [1] and additional stress in the Pd-matrix results in the 
solubility enhancements [2]. 

Solubility enhancements ratio (n /n) can be calculated from the equation [3]. 


n 

— = expl 
n 


( 


Oj^Vh 
RT 


1.1 


Where, 

Vh = Partial molar volume of hydrogen. 

T = Temperature of solubility measurement 
o = Tensile hydrostatic stress present in the Pd matrix. 

An and n are the atom fractions of hydrogen in equilibrium in the stress field. (ah m ) of 
the precipitate and in the absence of stress ( i.e. ah” - 0), respectively. 

In the present case above ratio is always greater than 1 .0 for Pd-AbCb composite 
because tensile stress are generated in the Pd matrix on cooling down. 
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Experimentally, we can find the solubility enhancement ratio. However, the procedure 
does not help in understanding how the strains and stresses are distributed within the 
composite. These can be easily studied by Finite Element Method (FEM). Here FEM has 
been applied to understand the stresses developed in Pd-Al 2 C>3 composites. The radial 
stress generated in the Pd matrix helps in calculating (n /n) ratio. 

1.2 Finite Element Method 

Finite element method (FEM) is a versatile method for the approximate solution 
of the differential equations. Its use in industry is enormous and it is said that without it 
many problems would remain unsol vable. FEM method is particularly applicable for the 
stress-strain relations as these relations involve differential equations. FEM has been 
applied for two types of models. Spherical precipitate is embedded in infini te spherical 
matrix in the first model and second model considers cylindrical precipitate in infini te 
cylindrical matrix. 

Hydrogen solubility enhancements in Pd-Al 2 0 3 composite only provides 
information about the maximum radial stress build in the Pd-matrix. It does not tell about 
the distribution of radial and tangential stress build in the composite. FEM model can 
provide the radial stress, tangential stress and tensile stress generated throughout the 
composite at different positions. 

1.3 Work Plan 

Elastic FEM formulation has been evaluated to study various stresses generated in 
Pd-Al 2 C>3 composite. Two types of models are evaluated through FEM formulation. 
Spherical precipitates in infinite spherical matrix are considered in first model whereas 
the second model deals with cylindrical precipitate in infinite cylindrical matrix. 

1.3.1 Plan For Spherical Precipitates In An Infinite Spherical Matrix 

Assumptions are made before applying this model. Respective force balance 
equations are established and stress equations are evaluated involving displacements and 
strains. The variation method was used for solving the FEM model. In order to simplify 
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the model, two noded element was used to get displacement at every nodal point. Strains 
and stresses were evaluated at every element to show how strains and stresses were 
generating in whole composite. Results were generated by varying the temperature by 
keeping the volume fraction of AI2O3 precipitates same, and the vice-versa. Finally, 
■ results are obtained for solubility enhancement ratio and Helmholtz free energy. 

1.3.2 Plan For Cylindrical Precipitate In An Infinite Cylindrical Matrix 

Firstly, basic assumptions were made before evolving this model. Respective 
differential equations were generated for stresses and strains. Variation method was 
applied to solve the differential equations. Axis-symmetric case was made to simplify the 
model, eight noded isoparametric serendipity isoparametric element was evaluated for 
generating displacements in co-ordinate axis. Stresses were transformed into polar co- 
ordinates from cartesian co-ordinates for calculating solubility ratio and Helmholtz free 
energy. 
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Chapter 2 


LITERATURE REVIEW 

2.1 Internal Oxidation And Hydrogen Solubilities In 
Composites 

If a Pd alloy containing a substitutional base element is exposed to an oxidizing 
environment at high temperatures, internal oxidation of the base element may occur if the 
solute is present in small quantities and it has been studied in a variety of noble metal- 
base metal alloy systems [4,5]. 

A1 particle internally oxidizes to AI 2 O 3 precipitates which results in formation of 
Pd- AI 2 O 3 composite at higher temperature. As the temperature is cooled from high 
temperature to room temperature, thermal residual stresses are generated at the composite 
due to mismatch of material properties of Pd and AI 2 O 3 precipitates in composite. 

Thermal residual stresses generation in composite is complex in nature. It is well 
known that Pd absorbs H 2 and this property of H 2 has been utilized to study stresses 
developed in the composite. Additional stress in the composite results in the solubility 
enhancements and given by equation 1.1. Hydrogen can be used as probe element to 
know the stress distribution in composite. Higher the solubility enhancements, higher will 
be stresses. These stresses developed depend upon the temperature of internal oxidation 
and precipitate volume fraction. 

2.1.1 Hydrogen Isotherms After Internal Oxidation 

Hydrogen solubility in Pd matrix helps in calculating hydrogen solubility 
isotherm in a-phase for alloy and pure Pd is shown in Figure 2.1. The Figure 2.1 shows 
that the low pressure solubilities decrease only slightly and the terminal hydrogen 
solubilities increase significantly with increase significantly with increasing mol fraction 
of Al, Xai. There is only a small difference between the dilute phase H 2 solubilities of Pd 
and the Pd-Al alloy because of the low Al content [1]. 
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The solubility of hydrogen at 323 K is enhanced in the dilute region of the 
internally oxidized Pd-Al alloys after complete oxidation at 1073 K, 1173 K and 1273 K 
as can seen in Fig. 2. 1 for Pd-Al alloys after they have been internally oxidized. 

In Fig. 2.1 the hydrogen content is expressed as the mol H-to-mol metal ratio 
where the latter is that of the initial, unoxidised alloy. If M is computed, instead, on the 
basis of mols of Pd matrix after oxidation, the calculated H/M values are very similar [1]. 
It can be seen that the hydrogen solubility enhancement decreases with increasing 
oxidation temperature. The solubility enhancement factor (n’/n) is defined as the ratio of 
hydrogen solubility at a fixed hydrogen chemical potential, i.e. fixed equilibrium 
hydrogen pressure, in the dilute phase after oxidation treatment (n’) to that in pure Pd (n). 
It has been shown [6] that the Pd matrix after internal oxidation of Pd-AI alloys is 
essentially pure, and therefore, if there were no influence of the ceramic particles on the 
solubility, the solubility enhancement, (n’/n), would be 1. 

The variation of the solubility enhancement factors alloys are shown as a function 
of equilibrium hydrogen pressures in Fig 2.2, which shows that the enhancement factor 
for complete oxidation at all of the temperatures is larger at the lower, as compared to the 
higher pressures where it becomes relatively constant. At all of the equilibrium hydrogen 
pressures, lower oxidation temperatures lead to greater solubility enhancements in the 
dilute phase isotherms at 323 K. 

2.1.2 Theoretical Analysis Of Hydrogen Solubility 

At the relatively high temperatures employed for the internal oxidation 983,1073, 
1173 and 1273 K, it seems reasonable to assume that the transformation stresses which 
result from the precipitation of alumina particles at the above temperatures are largely 
accomplished by the relaxation of the matrix. Stresses will develop, however, during 
cooling from the oxidation temperature. They will start to develop at a temperature below 
approximately half of the melting point (K) of Pd (913 K). The stresses develop because 
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Figure 2.1 -- Dilute phase hydrogen Isotherms for 
completely oxidized Pd 0 97 A1 0 03 measured after an initial 
isotherm measurement followed by evacuation at 323K: 
A, 1073K; 0, 1173K;pj, 1273K [ 1 ], 
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of the difference in the coefficients of thermal expansion of the palladium matrix (a m ) 
and the embedded alumina ceramic (etc) precipitates. Upon cooling, an imbedded ceramic 
particle contracts to smaller extent than the Pd matrix and hence the metallic matrix will 
be subjected to tensile and the ceramic precipitates to compressive stress. These states of 
stress resulting from cooling ceramic precipitates wi thin metallic matrices have been 
experimentally verified in analogous systems by X-ray diffraction, e.g. in Al-graphite [7] 
and Al-SiC particulate-composite [8] systems where an average tensile stress state was 
found in the matrix phase. 


The thermal stresses developed in the spherical alumina precipitates and matrix 
have been derived for composite systems by assuming elastic. Among the models which 
have been proposed, that of Kim et al, [2] seems to be the most realistic because it 
accounts for the volume fraction of the precipitates in a finite matrix, whereas the other 
models are based on a single ceramic particle in an infinite metallic matrix. According to 
this elastic model of Kim et al.[2] the tensile hydrostatic stress (ah™) developed in the 
matrix (superscript m), which can be shown to be independent of distance from the 
precipitate-matrix interface, is given by [1] 






2.1 


where Pi, the radial pressure at the metal-ceramic interface, is given by 



In the above equations, a™ is the thermal expansion coefficient of the matrix, etc is the 
thermal expansion coefficient of the ceramic precipitate, T is the difference between the 
temperature at which relaxation of the lattice ceases, which is taken as one-half of the 
melting point in Kelvin, and 323 K (the temperature at which the isotherms were 
measured), Em is the elastic modulus of Pd (123 GPa), E c the elastic modulus of alumina 
(380 GPa), v m the Poisson's ratio of the matrix (0.375), v c the Poisson's ratio of alumina 
(0.250), a is the radius of the precipitate and b is the average half distance between the 
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center of two neighbouring precipitates. For a dispersed spherical precipitate system, the 
relationship between a, b and cp (the precipitate volume fraction) is obtained from [9] 


b = a + 


a exp(8<p) 

6<p 1/3 



2.3 


The above expression is accurate at low volume fractions, i.e. when (p<=0.1. The 
tensile hydrostatic stress developed in the matrix is a function of the elastic properties of 
the matrix and precipitate, and the volume fraction of precipitate. It is not, surprisingly, a 
function of distance from the ceramic/matrix interface because the terms involving this 
distance cancel out when individual stress components are s umm ed to obtain the net 
hydrostatic stress[2]. The hydrostatic stresses in the matrix for different volume fractions 
of alumina, as calculated from equations (2.1) to (2.3), are tabulated in Table 2.1. 

If the tensile hydrostatic stress is much less than the modulus of Pd, which is true 
in the present case, then the interaction energy between the stress field and one mole of 
hydrogen atoms is given by [10] 


ri 

— = exp 
n 


RT 


2.4 


where Vh is the partial molar volume of hydrogen in Pd (Vh=1.77x10"° 6 m 3 /mol). The 
absorption of hydrogen in the stress field of the ceramic precipitate is given using the 
Boltzmann approximation as 


ri 

— = exp 
n 



2.5 


where n' and n are the atom fractions of hydrogen in equilibrium in the stress field (cO 
of the precipitate and in the absence of the stress field (i.e. Gh m =0), respectively. The 
above equation is integrated over the entire matrix from the ceramic-matrix interface (r 
=a) to half the distance of separation of the spherical ceramic precipitates (r =b). The b/a 
ratio can estimated from equation [2.3]. 


9 



<p 

a h/MPa 
(spherical) 

n'/n 

(spherical) 

<rh/MPa ’ 
(cylindrical) 

n'/n 

(cylindrical) 

0.02 

81 

1.06 

123 

1.09 

0.03 

148 

1.10 

127 

1.09 

0.04 

254 

1.18 

132 

1.00 

0.05 

298 

1 .22 

136 

1 .09 

0.06 

337 

1.25 

140 

1.10 

0.07 

372 

1.28 

145 

1.10 

0.08 

404 

1.31 

150 

1.10 

0.09 

433 

1.33 

154 

1.11 

0.10 

460 

1.35 

159 

1.11 


table 2.1 - Calculated hydrostatic stresses and 
solubility enhancement for spherical and cylindrical 
alumina precipitates in Pd for AT=590K and T=323K 


for n’/n. 


The enhanced concentration in the presence of the tensile stress field in the Pd- 
matrix due to the ceramic precipitates is then given by, 


fff expf-^- V 2 smQdrdQdfy 


n 

n 


\RT ) 


r r r ri 


= exp 


RT 


A 


2.6 


where E(a h m ) is given by equation (2.1) using equations (2.2) and (2.5); it is surprising 
that n'/n does not depend on r,0 and cp. This is because the hydrostatic stress developed in 
the matrix is only a function of a and b, and does not depend upon the distance from the 
ceramic-matrix interface. From the known elastic properties of the matrix and the 
precipitate and the volume fraction of the precipitates, the solubility enhancement in the 
matrix due tensile hydrostatic stress can be calculated using equation [2.6]. For the case 
of spherical alumina precipitates in Pd, the calculated solubility enhancements for 
different volume fractions of alumina are tabulated in Table 2.1. 


2.2 Residual Stresses 

Residual stresses are the system of stresses which can exist in a body when it is 
free from external forces. They are sometimes referred to as internal stresses, or locked- 
in stresses. Classically, residual stresses are categorized into three types according to the 
length scale over which they act [11]. Type-I stresses act over distances measured in 
millimeters and are often referred to as macrostresses. They can typically develop within 
the material following welding, machining, or surface finishing processes. Stresses that 
act over length scales that are typical of a microstructure are termed type-II stresses. The 
residual stresses that are developed in MMCs due to the interactions between the 
reinforcement and the matrix belong to type-II stresses. Finally, the stresses that exist at 
an atomic scale and vary over individual grains are described as type -III stresses. Type- 
II and III are often collectively referred to as microstresses. 
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2.2.1 Thermal Effects 

Thermal stresses in metal-matrix composites result during cooling either from 
temperature gradients within the specimen (leading to type-I stresses) or from a mismatch 
of coefficients of thermal expansion (CTE) between the constituents (introducing type-II 
stresses) [12]. The former effect can be controlled by slow cooling, but the latter is 
inherent and will always give interaction stresses unless the CTEs of the constituents are 
identical. Typically, the difference in CTEs in metal-matrix composites is large, for 
example Pd (cxm =1 1.5 x 1 0"° 6 /K) and AI 2 O 3 (ctc= 8.0 x 10"° 6 /K). Processing temperatures 
of metal-matrix composites are usually several hundred degrees above ambient. Thus, on 
cooling, the larger contraction of the matrix compared to the reinforcement leads to 
tensile and compressive residual stresses in the matrix and the reinforcement, 
respectively. In continuous fiber or aligned whisker composites, the anisotropic shape of 
the reinforcement induces larger axial than transverse residual stresses. 

Most of the early analytical models proposed to estimate the thermal residual 
stresses are based on Eshelby's equivalent inclusion method [13,14]. Using Eshelby's 
method, the stress fields in the region surrounding a reinforcement, which form as a result 
of a uniform temperature change, may be determined [15,16]. Eshelby's, method has also 
been used to predict the yield stress and the work hardening rate of metal matrix 
composites [17] as well as to estimate the overall thermal expansion characteristics of a 
heterogeneous solid containing a dispersion of aligned inclusions [18]. However, most of 
the published analytical studies are limited to simplified cases, such as those involving a 
uniform temperature change boundary condition, a single reinforcement in an infinite 
matrix, or the case of a perfectly elastic or elastic-perfectly plastic matrix. When the 
composite system involves non-homogenous deformation as a result of matrix plasticity 
and incorporates multiple reinforcements containing sharp comers, the analytical 
'formulation deduced from analytical methods can become intractable and the need for a 
numerical analysis arises. Numerical techniques, such as the finite element method 
(FEM), have been widely used by a number of investigators to predict the residual 
stresses and the constitutive response of the MMCs [19-24], 
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Thermal residual stresses have been studied in several metal-matrix systems: 
Al/SiC (Ho 1990 [19]; Ho and Saigal 1994 [21]; Levy and Papazian 1991 [22]; Shi and 
Arsenault 1991 [23]; Suresh and Nakamura 1993 [24], Ho and Saigal 1994 [25]; Suery et 
al. 1993 [26]), Ti alloys composites (Saigal et al. 1992 [27], Durodola and Derby 1994 
[28]; Bigelow 1992 [29]) and Pd alloy composites (Balasubramaniam et al. 1997 [1]; 
Eastman and Ruhle 1989 [30]), by a variety of techniques that include X-ray diffraction ( 
Tsai et al. 1992 [7]; Ledbetter and Austin 1987 [8] and Cox 1987 [31]) and neutron 
diffraction (Ledbetter and Austin 1987 [8]; Krawitz 1990 [32] and Majumdar et al. 1991 
[33]). 

Saigel et al. [27] measured thermal residual stresses developed during cooling of 
silicon carbide fiber-reinforced titanium aluminide matrix composite using neutron 
diffraction and compared with those obtained by finite element analysis. Figure 2.3 
illustrates the estimated thermally induced residual stresses in the titanium aluminide 
matrix as the composite is cooled from strain- free temperature of 806°C. 

In modeling, the process of generation of thermal residual stresses in cast SiC 
fiber-reinforced aluminum alloy A357.0 matrix composites, Ho and Saigal [19] employed 
thermo-elastoplastic finite element analysis. Figure 2.4 shows the von Mises’ equivalent 
stresses in the matrix around fiber at the fiber/matrix interface as the temperature 
decreases from 620 to 27°C [19]. 

Suery calculated the thermal residual stresses in 15 vol.-% SiC particulate 
reinforced Al 1100-0 composite [26]. In their study, it was found that the reheating of 
the composites after quenching to room temperature changes the sign of the radial stress, 
which become positive (Figure 2.5). Figure 2.5 shows the variation of the radial stress 
with temperatures (from 660°C down to 20°C and up to 660°C) at the particle/matrix 
interface. The combined effects of thermal residual stresses and fiber spatial distribution 
on the deformation of a 6061 aluminum alloy containing a fixed concentration 
unidirectional boron fibers have analyzed by Suresh and Nakamura [24] using detailed 
finite element methods. The model involves an elastic-plastic matrix, elastic fibers, and 
mechanically bonded interface. The results indicate that both the fiber packing and 
thermal residual stresses can have a significant effect on the stress-strain 
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Figure 2.3 - Estimated residual stresses induced in the 
Ti-14Al-21Nb matrix reinforced with 35% SCS-6 fibers 
cool down from the strain-free temperature (806° C) to 
room temperature [27]. 
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Figure 2.4 - Equivalent stress in the A357.0 matrix 

around SiC fiber at the fiber/matrix interface ( air 
cooling: temperature of air is 27° C) [19]. 





Figure 2.5 - Radial stress at the interface as a function 
of temperature during cooling to room temperature at 
K/s and subsequent heating at various rates [26]. 


<£> 




characteristics of the composite. The thermal residual stresses cause pronounced matrix 
yielding which also influences the apparent overall stiffness of the composite during the 
initial stages of subsequent far- field loading along the axial and the transverse [21]. 
Thermal residual stresses developed during casting of SiC/aluminum particulate- 
reinforced composites were investigated by S.Ho and A.Saigal as a function of cooling 
rate and volume fraction of particles using thermo-elastoplastic finite element method 
[ 22 ]. 

It was found that the matrix undergoes significant plastic deformation during cool 
down and has higher residual stress distribution as the cooling rate increases. The 
presence of thermally induced residual stresses tends to decrease the apparent modulus of 
elasticity and increase the yield strength of the composite compared to those without 
residual stresses [22]. 

2.3 Brief Introduction Of Finite Element Method (FEM ) And 
Its Application 

Finite Element Method (FEM) is a powerful tool for the approximate solution of 
differential equations governing diverse physical phenomena . Its use in industry & 
research is extensive, indeed it could be said that without it many problems would be 
incapable of solution [35]. This method is normally applied when conventional method 
fails. Finite Element methods originated in the field of structural analysis, and many of 
the early applications were in civil and aeronautical engineering. Methods have also been 
developed for all types of two- and three-dimensional solid mechanics and recently have 
been extended into fluid mechanics, composites and heat transfer fields [34]. 

All finite element methods involve dividing the physical system, such as 
engineering componenfcor structure, into small sub-regions or elements. Each element is 
an essentially simple unit, the behaviour of which can be readily analysed. The 
complexities of the overall system are accommodated by using large number of elements. 
In the case of a pin-jointed structure, for example, the choice of elements is clear; each 



member of the structure is treated as an element, subject only to a tensile and 
compressive force. From simple analyses of the stiffness of the individual elements, it is 
possible to establish the behaviour of the overall structure in terms of, say, the 
displacement responses to externally applied loads. These loads are often applied at the 
joints of the structure, and it is convenient to consider the displacements at the same 
points, which are usually referred to as nodes or nodal points [35J. 

Frameworks in which relative rotations of the elements at the joints are restricted 
can be treated in a similar fashion except that, in addition to forces and linear 
displacements, bending moments and rotations must also be considered. When attention 
is extended to structure involving plates and shells, however, the choice of elements is 
less obvious; the continuous surfaces must be sub-divided. Nevertheless, once the 
elements have been chosen, the behaviour of each can be analysed in terms of 
displacement response to applied forces and moments, the variables being referred to the 
nodal points at which the elements are effectively connected together. Similar remarks 
concerning the choice of elements and the significance of nodal points can also be made 
in connection with other two- and three-dimensional solid continuum problems[34]. 

At first sight, some finite element methods may appear to be fundamentally 
different from analytical and finite difference methods in that they do not use the 
governing differential equations explicitly. They are, however, based on the physical 
principles of equilibrium, compatibility and constitutive relationships used to derive the 
differential equations. 

2.4 Helmholtz Free Energy 

There are four types of thermodynamic potentials but Gibbs energy (free 
enthalpy) and Helmholtz free energy are the most important. The free energy is the 
portion of the internal energy available for doing work at constant temperature. The 
enthalpy h as defined here is the portion of the internal energy that can be released as 
heat when thermodynamic tensions are held constant. The internal energy density u is a 
potential for the thermodynamic in an isentropic process (s = const.), while the 



Helmholtz free energy density calculates potential for temperature and equation is given 

in terms of strain and temperature:- 

¥ = jj£ y£ ji +~^ 2 -(3A + 2ju)ATa.e 

2.7 

where e = £„. + s ee + for polar coordinates. 

2.8 

where e = + e a for cartesian coordinates. 

2.9 
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Chapter 3 


SPHERICAL PRECIPITATES 

In the first model that is considered for evaluating stresses and strains in the Pd- 
AI 2 O 3 composite, the spherical AI 2 O 3 precipitate is embedded in infini te spherical Pd 
matrix. It is cooled from high temperature to room temperature, pure elastic case is 
considered for simplification. This method was developed by Lee [2] and is popularly 
known as Lee model. It is necessary to consider assumptions before evaluating any 
mathematical model. Assumptions are basically the special cases which help to simply 
the mathematical model. 

A numbers of assumptions are made to evaluate this model [2]. The reinforcement 
AJ 2 O 3 ceramic precipitate is considered as a sphere in an infini te matrix of metal-matrix 
composite (MMC). The matrix behaviour is elastic in nature . The stress-strain behaviour 
is independent of strain-rate and stress oriented. The temperature in MMC is uniform at 
all times in all directions. Matrix containing spherical AI 2 O 3 precipitate is shown at 
Figure 3.1. For mathematical simplification, we assume that all other ceramic precipitates 
are of spherical shape with radius a, and further, that a given particle is surrounded by 
uniformly distributed neighbouring particles [ 2 ]. 

This assumption allows us to use spherical coordinates (r, 0 and cp) whose origin 
is at the center of a ceramic particle with radial symmetry as shown in Figure 3.1. 

3.1 Formulation Of Basic Equations 

Consider a thick hollow sphere whose internal surface of radius a contains 
ceramic particle and outer surface contains matrix as shown in figure 3.1. Because of the 
radial symmetry approximation, the tangential displacements as well as the shear stresses 
are all zero [2]. The radial displacements u is a function of radial distance r. Further, the 
force balance in r direction in the absence of body forces leads to the following 
equilibrium equations. 
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Figure 3.1 - A schematic diagram for an idealized 
ceramic configuration in Infinite spherical matrix. 



3.1 


for , 2 (0~r-^) _ q 
^ r 

where cr r is radial stress and <j g is tangential stress 

The radial and tangential strains are defined as follows 

s= — (radial strain) 3.2 

dr 

e g = — (tangential strain) 3.3 


Assumption is made that the stress and strain in third direction (<p) is equal to 
stress and strain in tangential (0) direction. So the stress and strain in 3 rd direction is 
given by 

<j 6 — and £ q £ ^ 

Under these conditions the necessary stress-strain temperature relations are 


s r = -£-(cr r - 2ua g )+ a AT 
Jb 

S e = 4r[- V(J r + (l - u)crg]+aAT 


3.4 

3.5 


where E = Young’s modulus of particle. 
o = Poisson’s ratio of particle. 
a = Coefficient of thermal expansion of particle. 

A T = Temperature difference between high temperature and room temperature. 
Equations 3.4 and 3.5 can be inverted to give the radial and tangential stress in 
terms of the strains and the temperature. 


E, 


CTg'] 2 j (l + uXl _ 2t>) 


1-t; 2 

oV2 1 


_ KaKT_ 
s a S\ (l + ^Xl-2u) 


l-t» 0 V 2 
v4i 1 
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[flL= (1 + 4 - 20 ) 

where [D] e is known as the elastic property matrix which is a symmetric matrix. 

Where [D]xh = [D] e a .AT is known as the thermal property matrix which is a 
symmetric matrix. 

3.2 Variational FEM Method 

An alternative, although entirely equivalent, approach to that of solving the 
governing differential equation is the variational method which seeks to solve a problem 
by finding the conditions for a stationary value (either a maximum or a minimum) of 
some quantity with respect to the unknowns of the problem [23]. This quantity, which is 
often referred to as a functional, generally involves an appropriate integral of the 
unknowns or derivatives of the unknowns over the region of interest, which in 
engineering elasticity problems is the component or structure concerned. Variational 
formulations for elasticity problems can be developed from either the governing 
differential equations or an equivalent physical principle, such as an energy minimization 
[23]. Here the first approach is followed. The equilibrium equation (3.1) leads to 


1-u v-Jl 
uV2 1 


| Ua r 

+ 2K - 


.5u.4nr 2 .dr = 0 

o J l dr 

r 

J 


) (do r 

+ 2ct '- 

2<V 


JyuAitr 1 .dr = 0 

o l dr 

r 

r , 



After integration by parts, the matrix form of this equation becomes 





' du 

dr \dr = 0 

V2 u 


3.8 
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Radial Distance 



Figure 3.2 — The discretization of the domain for the 
case where AI 2 O 3 precipitate is inside the composite. 
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Substitution of equations (3.6), (3.2) and (3.3) and use of properties of variational 
operator leads to the following equations 

We get 51 = 0 3.9 


where I = ] ]-r l \sf [d \ {*}- r 2 [dJ, {f} 
o L z 


.dr 


3.10 


where I is the variational functional of the problem. 


3.2.1 Shape Functions 

The variation of the dependent variable over each element is often expressed in 
terms of the variable at certain points called as nodes and the known functions are called 
shape functions. The discretization of the domain of the composites shown on Figure 3.2. 
In the present case a simple approximation consisting of 2 nodes and linear shape 
functions is used : 

e 

1 2 

Ui U e 2 

simple element for evaluations of shape functions 
u = N'u{ +N 2 u‘= {w} e 3.11 

where u is the displacement term with respect of shape functions (N} e and elemental 
displacement {u} e . 

where N? = — — — and N 2 =— — — 3.12 

r \ ~ r 2 r 2 T \ 

The derivatives of shape functions with respect to r are given by 


dN[ 1__ 

dr r“ - r 2 


and 


dr 


r'~, - r, e 


3.13 


25 



Radial strain in terms of derivatives of shape functions and displacements is given by 

3.14 


du \dNY' Xa 

*■ V*J M 


Strain vector terms of displacements and shape functions from equations (3.11) and 
(3.14) is given by 


■W- 


■W'W 


3.15 


\ s rr 1 idu/dr 1 dNjdr dN 2 /dr 
[V2 s ffe J~ |V2 u/r\ ~ (_V2 N ] /r JlN 2 /r\\u 2 

Here, [6] e is the elemental displacement vector containing the radial displacement at each 

node, substitution of equation 3.15 in the expression for I leads to 

/ 


3 * 2 


M = [D\{B}‘ .dr 


3.16 

3.17 


where [k] e is the elemental coefficient matrix and is obtained from the 1 st integral of 
equation 3.16. 

•i 


\ft = ymDXMdr 


3.18 


where {f} c is the elemental right side vector containing the thermal residual stress term. 

The elemental displacement vector {8} e is related to the global displacement 
vector { A } by the following relations 

^} e =M tf {A} 3.19 

where the matrix [A] e depends on the elemental (i.e. local node) and global node 
numbering systems. Substitution of equations 3.19 and 3.14 gives 


em\ ^ 

[k] = ^ [a] T [/t] e [a]" where [K] is the global coefficient matrix. 


3.20 


3.21 


tf-1 
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where {Frh}is the global right side vector. 


3.22 


{f n }=tww' 

Equation 3.20 leads to 

/^{aFMaMfKa} 

Extremization of I leads to 

M{A}={F ra } 


Solving of the above equations gives displacements at all the nodes. 
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Data for FEM Formulation 


Composite consists of two materials, Palladium (Pd) is the matrix and alumina 
(AI 2 O 3 ) is the ceramic precipitate & their material properties is given below in table 3.1. 


Criteria 

Material Properties of Pd- 

matrix 

Material Properties of 

AI 2 O 3 Ceramic particle 

Poisson’s Ratio (v) 

0.375 

0.25 

Young Modulus (E) 

123 * 10 y Pascal (N/m^) 

380 * 10 y Pascal(N/m 2 ) 

Coefficient of thermal 

Expansion(CTEs) 

1 1.5 * 10^/K 

8.0 * lO^/K 

s 


Table 3.1 - Material properties of composite. 


Calculation of [D] and {T} matrices 

For these materials, the matrices appearing in equation (3.6) are as follows : 

For Pd-matrix 


[4» = 


4 E„ 


11 


5 3V2 
3V2 


8 


[DL =11.5*10^ AT [DU 


For AI 2 O 3 Ceramic Precipitates 



Ey* [l2 4V2 

10 4V2 10 


[dIt, = 8.0*10' 6 AT [DJppj 
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3.3 Estimation Of Displacements 

Displacement of node can be evaluated by FEM formulation is given by solving 
equation 3.10 and theoretical displacements is given by 

u = arAT 3.24 

Validation of results is made by using only one type of material, it may be total 
AI 2 O 3 precipitate or may be total Pd-matrix. Exactly same type of results are obtained for 
both FEM modeling and theoretical equations and it is shown on table 3.1. Two kinds of 
cases are considered in analyze the FEM formulation result and shown in Figure 3.3. In 
the first case called 10 case where AI 2 O 3 ceramic precipitate is inside and Pd matrix is 
outside. In the second case called 01 case where ceramic precipitate is outside and Pd 
matrix is inside. 

Thermal residual stresses occur due to temperature difference and the temperature 
difference is kept at 1000K at both the above cases to obtain the various results. Figure 

3.4 shows the displacement for the composite as a function of radial distance. The Pd 
matrix has higher displacement because its coefficient of thermal expansion (a m = 1 1.5 x 
lO'VK) is higher than that of AI 2 O 3 ceramic precipitate (a c = 8.0 x 1 0 ' 6 /K). The change 
of slope at the interface occur due to the presence of AI 2 O 3 and Pd matrix at the interface. 
10 and 01 case are taken from the Figure 3.3. 

3.4 Estimation Of Strains 

Radial and tangential strains are calculated from equation 3.15 . 

Where 8 rr is the strain at each element at radial direction and see is the strain at each 
element at the tangential direction. 

Figure 3.5 shows the variation of radial strain for AT = 1000K for 10 and 01 
cases. The figure shows the radial strain is higher for Pd matrix (otm = 1 1 .5 x lO^/K) then 
that of AI 2 O 3 precipitate (etc = 8.0 x 10^/K). No temperature changes occur within 
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Radial Dist. FEM 

Theoritical 

FEM Theoritical 

In Mts. 

Calculation 

Calculation 

Calculation Calculation 

Total Pd Matrix 


Total 

Ceramic Particle 

0 

0 

0 

0 

0 

0.001 

1.15E-05 

1.15E-05 

8.00E-06 

8.00E-06 

0.002 ' 

2.30E-05 

2.30E-05 

1.60E-05 

1.60E-05 

0.003 

3.45E-05 

3.45E-05 

2.40E-05 

2.40E-05 

0.004 

4.60E-05 

4.60E-05 

3.20E-05 

3.20E-05 

0.005 

5.75E-05 

5.75E-05 

4.00E-05 

4.00E-05 

0.006 

6.90E-05 

6.90E-05 

4.80E-05 

4.80E-05 

0.007 

8.05E-05 

8.05E-05 

5.60E-05 

5.60E-05 

0.008 

9.20E-05 

9.20E-05 

6.40E-05 

6.40E-05 

0.009 

0.0001035 

0.0001035 

7.20 E-0 6 

7.20E-05 

0.01 

0.000115 

0.000115 

8.00E-05 

8.00E-05 

Table 

3.2 ~ Validation of results with theoretical and 

FEM results. 






Radial Distance 


10 Case 


01 Case 



Figure 3.3.a 

Ceramic ppt. is inside and 
Pd is outside at AT=1000K 
(ceramic volume 
fraction = 0.125) 



Figure 3.3.b 

Ceramic ppt. is outside and 
Pd is inside at AT=1000K 
(ceramic volume 
fraction = 0.875) 
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10 Case 01 Case 



Figure 3.3.a 

Ceramic ppt. is inside and 
Pd is outside at AT=1000K 
(ceramic volume 
fraction = 0.125) 


Figure 3.3.b 

Ceramic ppt. is outside and 
Pd is inside at AT=1000K 
(ceramic volume 
fraction = 0.875) 
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Displacements in metres 



Figure 3.4 -- Variation of displacements as a function of 
radial distance at AT = 1000K for IO and OI cases. 
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Radial strain e 



Figure 3.5 ~ Variation of Radial strain (e rr ) as a 

function of radial distance at AT — 1000K for IO and 
01 cases. 
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element during FEM modeling. Radial strain at initial part upto interface is infl uenced by 
individual part and temperature plays no part in influencing radial str ain . That’s why 
radial strains remain parallel to x-axis for initial part. At the junction of interface, slope of 
radial strains change drastically due to the mismatch of material properties. 

Figure 3.6 shows the variation of tangential strain as a function of radial distance 
AT = 1000K. Pd matrix (a m = 11.5 x lO^/K) has higher strain than AI 2 O 3 precipitate (etc 
= 8.0 x 10 6 /K) due to higher coefficient thermal expansion of Pd matrix. Slope of strain 
changes drastically at the interface due to mismatch of material properties of AI 2 O 3 
precipitate and Pd matrix at interface. 

3.5 Estimation Of Stresses 

Radial and tangential stresses are calculated from the strain in each element using 
equation 3.6. The stresses are depending on strains, temperatures and material property of 
the particles. Radial stress variation as a functions of radial distance at AT = 1000K is 
shown at Figure 3.7. The ceramic volume fraction of alumina precipitate is 0.875 for 01 
where as for 10 case is 0.125. The Young's modulus of alumina (E = 380 Gpa) is more 
than that of Pd matrix ( E = 123Gpa), stress is greatly influence by alumina than Pd 
matrix and has high negative value for 01 case than that of 10 case. The 10 case radial 
shown value has less negative value due to less presence of alumina particle in the 
composite. 

The stresses remain constant for initial part as the particle is solely influenced by 
individual particle change and no temperature changes occur within the composite during 
FEM formulation. Figure 3.8 shows the variation of tangential stress as a function of 
radial distance at AT = 1000K. The interesting thing about this graph that stresses in 
alumina is always positive whereas stress in Pd is always negative. The reason 
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Tangential Strain e 



Figure 3.6 — Variation of tangential strain (see) as a 
function of radial distance at AT = 1000K for IO and 
OI cases. 
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Figure 3.7 -- Variation of radial stress (a rr ) as a 

function of radial distance at AT = 1000K for IO and 
OI cases. 
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Figure 3.8 -- Variation of tangential stress (dee) as a 
function of radial distance at AT = 1000K for IO and 
OI cases. 
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behind this, may be that alumina tangential stress must have tensile stress as it has high 
Young’s modulus value. To compensate the tensile value of alumina precipitate, Pd 
matrix must have compressive stress due to less Young's modulus value. 


3.6 Estimation Of Solubility Enhancement ratio 

Hydrogen is only observed by Pd matrix and not by the alumina ceramic 
precipitate, This property of Pd is used to calculate thermal residual stress. According to 
Kim et al. model [2], the tensile hydrostatic stress ( Ch m ) developed in the matrix which 
can be shown to be independent of distance from the precipitate -matrix interface, is 
given by equation 2.1. The tensile hydrostatic stress ( Oh" 1 ) curve which is estimated from 
the radial stresses is shown in Figure 3.9, the radial stresses being considered here is not 
at metal-ceramic interface but over the whole composite. 01 case has more negative 
radial stress , so it has more positive tensile hydrostatic stress curve. 

The solubility enhancement (n /n) ration can be calculated after having calculation 
of tensile hydrostatic stress (ah m ) and is given by equation (2.5). The solubility 
enhancement value changes with the changing value of ceramic volume fraction at 
particular temperature difference. These solubility enhancement values are plotted in 
Figure 3.10 and Figure 3.11 as a function of alumina volume fraction for different 
temperature difference for 10 and 01 respectively. Figure 3.10 shows ratio as one for zero 
alumina ceramic fraction, ratio goes to maximum with the increase of fraction. Then with 
the increase of ceramic fraction the ratio decreases slowly and becomes one at 100% 
ceramic particle. The Figure 3.10 clearly shows the higher temperature difference as 
higher solubility ratio. Figure 3.11 show solubility enhancement curve for 01 case. It 
starts with value 1 and with the increase of ceramic fraction, ratio goes on increasing. 
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Tensile Hyd 
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Radial distance (metres) 


Figure 3.9 -- Variation of Tensile Hydostatic stress 
(a h m ) of as a function of radial distance at AT = 1000K 
for IO and OI cases. 


38 



Al 2 0 3 volume fraction 


Figure 3.10 — Variation of solubility enhancement 
ratio as a function of AI 2 O 3 precipitate volume fraction 
for different temperature difference values for IO 
case. 
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Solubility enhan 



Figure 3.11 -- Variation of solubility enhancement 

ratio as a function of A1 2 0 3 precipitate volume fraction 
for constant temperature difference for 01 case. 
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3.7 Estimation of Helmohltz Free Energy 


Helmohltz free energy density estimated potential of the particle for the 
temperature changes. Helmohltz free energy per unit volume in polar co-ordinates can be 
calculated from equation 2.7 as . 


V= 


2(1 + v) 


+ 2e 


22 


96 


M 


Eo 


2 (1 + u)(l - 2o) 


i £ rr + ) 


( 1-2 u )^ +2s ^ 


3.25 


The Helmohltz free energy over the whole volume is given by 

b 

\|/= J y/Anr 2 dr 3.26 

0 

Helmohltz free energy (HFE) depends on material properties, temperature 
changes and strains. Figure 3.12 shows HFE curve as a function of radial distance at AT = 
1000K. In both cases, energy for Pd matrix is more than alumina precipitate energy 
obtained by Pd matrix 10 case is very much as compared to other energies. At particular 
ceramic volume fraction, summation of HFE is done for alumina precipitate and also Pd 
matrix. Total HFE is obtains just by adding total HFE for alumina and Pd matrix. 

HFE values changes with the changing value of alumina volume fraction for 
particular temperature difference. These values of HFE are plotted as a function of 
ceramic volume fraction for 10 and 01 case at Figure 3.13 and 3.14 for different 
temperature difference. The HFE for alumina precipitate is very less as compared to Pd 
matrix. The Figures 3.13 and 3.14 clearly shows that higher free energies are obtained for 
higher temperature difference. The alumina precipitate HFE is zero at initial part where 
Pd matrix HFE has maximum value. Alumina precipitate HFE increases with the increase 
of ceramic volume fraction and goes to maximum at 100% ceramic volume fraction. Just 
the opposite case occur for Pd matrix, it has lowest value at 100% ceramic fraction. The 
HFE is mainly influenced by the Pd matrix as the percent of Pd decreases with the 
increases of al um ina precipitate, total Helmohltz free energy shows decreasing value with 
increasing ceramic fraction. 
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Helmholtz free energy (Joules) 


500 



Figure 3.12 » Variation of Helmholtz free energy as a 
function of radial distance at AT = 1000K for IO and 
OI cases. 
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Figure 3.13 -- Variation of Helmholtz free energy as 
a function of AI 2 O 3 precipitate volume fraction for 
different temperature difference values for IO case. 
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Helmholtz free energy (Joules) 


At AT= 590K 


At AT= 1 000 K 
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Figure 3.14 — Variation of Helmholtz free energy as a 
functions of AI 2 O 3 precipitate volume fraction for 
constant temperature difference for OI case. 


44 



3.8 Limitations Of The Model 


The spherical i n f in ite matrix considers only spherical precipitate and is based on a 
number of assumptions. These assumptions are basically the short coming of model. The 
major drawbacks of these model that no temperature change occur within the composite 
and model based on elastic case. Model gives the results for r and 0 direction and fail to 
give results for third component <p, rather it assumes 0 = cp due to radial symmetry. 

New model for cylindrical precipitate has been evaluates in the next chapter. It 
considers cylindrical precipitate in infinite cylindrical matrix (plane strain condition). The 
model considers temperature change within the composite. 
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Chapter 4 

CYLINDRICAL PRECIPITATES (CASE - 1) 


AI2O3 precipitates are embedded in infini tely long cylindrical Pd matrix. It is 
cooled from high temperature to room temperature. Numbers of assumptions are made to 
evaluate this model. Basic assumptions of this model is that strain in z direction is 
neglected, so it is case of plain strain condition. The distance taken in x and y direction is 
same. The matrix has pure elastic behaviour. The stress-strain behaviour is independent 
of strain-rate and stress oriented. The temperature in MMC is vary all the directions in 
within the composites. The temperature near the origin changes drastically. For 
mathematical simplification, we assume that a given particle is surrounded by uniformly 
distributed neighbouring particles. This assumption allows us to use spherical coordinates 
whose origin is at the center of a ceramic particle with radial symmetry as shown in 
Figure 4.1. Each composite quarter element is divided into eight (8) divisions as shown 
in Figure 4.1, 1 st division starts from x-axis and 8 th division ends on y-axis. 

4.1 Formulation Of Basic Equations 

Consider an infinite cylindrical matrix whose internal surface contains ceramic 
particle and outer surface contains Pd matrix. The domain consists of four quarter parts 
and only first quarter is considered to study the FEM formulation as shown in Figure 4.1 . 
Figure 4. 1 shows the discretization of the domain for the case where AI2O3 precipitate is 
inside and Pd matrix is outside. 

Any body subjected to external forces (body forces, surface forces) is described 
by the equations: 

(i) Variables 

Displacements: u, v, w which are functions of (x, y, z). 

Strain components: £>«, £yy, &zz, Sxy, £yz and Szx- 

Stress components: CTxx, ayy, CTxy, errand 

Where x, y and z are the direction of cartesian coordinates. 
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(ii) In the present case, basic assumption is made that strain in third direction 
(z ) is zero, so it becomes a plain strain case. Stress-strain relations for 
plain strain case is given below: 

Strain in z direction is zero. 

e y*=£ a =e a =0 

= 0 and <J yz =2^e y2 =0 


= 4e a + s yy )+ 2 fjs a - fa + 2/u)ocT 4. 1 

+ S yy )+ 2 ^yy ~ fa + 4.2 

= a(s„ + £yy )+ 2fis B - fa + 2faaT 4.3 

°xy ~ 2 H S xy 4.4 


where A = — and a = G = — r 

(l + v)(l - 2v 2(l + v) 

v = Poisson's Ratio of material. 

E = Youngs Modulus of material 
T = Temperature of the composite of material. 

H = G = Shear Modulus of material 
a = cofficient of thermal expansion of material 

4.5 

4.6 
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The strain equations are given by 

The strain equations are given by:- 

e„ - — = strain in x direction 
dx 

dv ...... 

£„ = — = strain in y direction 
w dy 



Distance in y—axis 


A 


7th Div. 


th Div, 


O represents element number 
| | represents nude number 
5th Div. =A^ Ceramic precipitate 

I — | = fd matrix 


4th Div. 



3rd Div. 


2nd Div. 


1st Division. 


0 E 


S 1 ® 0 


Distance in x-axis 


Figure 4.1 - Discretization of the quarter domain using 
2-D elements. 


4.7 


s = — (— + — ) = shear strain in x-y direction 
1 2 oy ox 

Equations (4.1) to (4.4) can be inverted to gives the stresses in terms of strain (from 
equations 4.5 to 4.7) and temperature in matrix form: 




► = 

X + 2/i 

X 

X 

X + 2/u 

O' 

0 


>-{3X + 2{i)ccT< 

V 

1 




0 

0 




0 

V J 


4.8 



X + 2ju 

X 

0 


X 

X + 2 /j. 
0 


0 

0 


where [D] e is known as elastic property matrix [35] . 


4.9 


It is clear that [D] e is a symmetric matrix. 

rn 


[■Dk =(3A + 2/^)a] 


1 

0 


> where [D]xh thermal elastic property matrix. 


4.10 


4.2 Variational FEM Method 

Variational approach has been discussed in section 3.2 for evaluating infinite 
spherical matrix. Variational formulations for elastic problems is developed either from 
governing differential equations or an equivalent physical principle. Second approach is 
followed to formulate the FEM problem. The strain energy in absence of body force is 
given by: - 

Strian energy = I = + <r„s„ + 2 cr^s^dA 4. 

n ^ 
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Strain energy in matrix form of stress and strain is given by 

4.12 

p ^ 

Putting equations 4.8 in 4.12 gives 

4.13 

The variation of the dependent variable over each element is often expressed in terms of 
the unknown values of variable at certain point called a node and the known function is 
called the shape functions. Here, 8-noded isoparametric serendipity element is used to 
evaluate shape functions and mapping part element is shown on Figure 4.2. 



Mapping Approximation 


Figure 4.2 ~ 8-noded isoparametric serendipity 
element. 




The displacements and temperature in terms of shape functions are given by: 


u = N[u{ + N e 2 u' 2 + 

+aw=WM* 

4.14 

V = N‘v’ + N[v\ + 

+jv;w=WW’ 

4.15 

T = N'T' + N'T' + 

+ Ar 8 T 8 c -{ATO' 

4.16 


Where u and v are the displacements in x and y directions only and T is the temperature 
term within the composite. Temperature in the composite is given by 


Temperature = Room temperature + temperature difference 


' b ' 
log— 
r_ 

, b 
log— 

V a) 


where b =radius of precipitate in composite, a = radius at 1 st node and r = radius. 


The individual shape functions is given by: 

Nl ='j(l + fXl-’7X- 1 + #-vX 

tf J *=l(l + fXi + ’7X-l+£ + 7). 

N‘=U i-fXi+’/X-W+i) . 


X*-v)- 4.17 


Putting the values of derivatives of shape functions and element displacements, we obtain 
strain matrix from equations 4.5 and 4.7 


.tNTRAL UHWT 

1 1. T«, KWi 
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du 











4.2.1 Calculation Of (B} e Matrix 

In order to calculate coefficients of {5}' , the derivatives of the shape functionwith 
respect to x and y are required, while these functions are given in terms of the intrinsic 
coordinates [35]: - 

dN‘ j dx dy 1 \dN' 

. L ^ 419 

dN[ dx_ dy_ QNJ_ 4 ’ iy 

drj j L 5; 7 3«Jl dy . 

'fflf 

~ [/} * > where [J] is the jacobian matrix. 

. ^ - 

Inverting equation (4.18), the required shape functions derivatives can be obtained 

' dN[ ) \ dN' 

dx _ r rl-l I 
' dN; bn: 

„ dy J dr] 

[J]' 1 is impossible to solve as we cannot express shape functions in x and y rather 
we can express shape functions {N} e in terms of £ and rj . 

We can easily evaluate [J] 





dx dy_ 

d£ dt 

dx dy_ 
drj drj 



'Em;-*; 

1 



/-i 


drj 


E dN ‘y, 

_M 

t dN ;y 

i-i 

drj 


dN' dN* 
d£ d£ 
dN[ dN' 
dr] dr] 


dNl 

dNl 



y 

y 2 

y $ 


Putting equation (4.19) into eq. (4.18) gives {B} e . 

{ 5 } e is the element displacement vector containing radial displacement at each 
node, substituting equation (4.19) for I leads to 


= jf WW - 

A 2 


dA 


1=1 


/« 1 L 


i VfmY-bfV»Y 


•where [k]* = ^{b)* T \ p\{B} e dA = elemental coefficient matrix 

n 

{ f n } = |{7’ } e [Dl {5}' dA = elemental right side vector 

n 


4.20 

4.21 

4.22 


4.2.2 Gaussian Quadrature 

With higher-order elements, in which not only displacements but also str ains and 
stresses vary over each element, the integration involved in the variational formulations 
generally cannot be evaluated analytically [35]. If there is no restriction on the positions 
at which the functions can be evaluated, then there is a class of numerical integration 
procedures refereed to as Gaussian quadrature which to be preferred. 

For integration of a functions, f( £, 77 ), of the independent variables over the two- 
dimensional region defined by - 1 < £ < +1,-1 < 7 < +1 . The integration is given by 


53 



4.23 


+ 1+1 n. n. _ 

-l-l M M 

where the , w, and tj- , w, are the pairs of abscissae and weight factors 

established for functions of single variables [35]. 

Equation (4.21) is given by: - 


{f} n= W T [ D l W det I J 

i-i >i L 



n-ij 


Where det |J| is the determinant of jacobian matrix. 


Equation (4.22) is given by 

{/V = i] f] w, w jf[{ T Y IPb {b}' det I J 

i-i j = i i-ij 

We have used four noded gauss values to evaluate integration and their values are 
§!=rn = -0.86113, wi = 0.347855. 

& = ri 2 = -0.33998, w 2 = 0.652145. 

& = ri 3 = 0.33998 , w 3 = 0.652145. 

= ri 4 = 0.861 13 , w 4 = 0.347855. 

The elemental displacement vector {<?} e is related to global displacement vector 
{a} by the following equation: 

W=WM 4.24 

where the matrix [A] e depends on the elemental (i.e. load) and global node 
necessary systems. Substitution of equation 4.22 in 4.19, we get 

^WMaMaFIf] 4.25 

w-st* n* ]* [Aj where [K] is the global body force vector. 

rt C 

where [F] is global right side vector. 

e=\ 


54 



4.26 


Extremization of I leads to (i.e. — = 0 ) 

dA 


Solving the above equation gives displacement at all the nodes in X and Y 
directions. 


Calculation of [D] e and [D]xh matrices 

Composite consists of Pd and A1 2 0 3 materials and there [D] e and [D]ih matrices 
are different for different material. 

For Pd-matrix 


[£>], = 4E -‘ 


11 


5 3V2 
3V2 8 


11.5*10-“ 


For AI2O3 Ceramic Precipitates 


[z>L=— 

1 * 10 


12 4V2 
4V2 10 


[Z)L=2^.8. 0*10 


cc 



4.3 Estimation of Displacements 

Displacements at all nodes is the x and y direction can be solved by solving 
equation 4.26. Here, u means displacement in x direction and v means displacement in y 
direction. The composite is made to cool from high temperature to room temperature, in 
discussing the model AT = 800K means the difference between the temperature inside 
the composite to the outside the temperature is 800K. The temperature profile is shown at 
Figure 4.3 for AT = 800K. Two types of cases is considered for analyzing FEM model 
and these are shown in figure 4.4. The first case will be called 10 case where AI 2 O 3 
precipitate is inside the system and Pd matrix is outside and in the second 01, the AI 2 O 3 
precipitate is outside and Pd matrix is inside. 

Figure 4.5. a and Figure 4.5.b show the variation of displacement of nodal points 
in x-axis line for 10 and 01 case at AT = 800K as a function of x-axis. In 10 case, 
displacements in x direction (u) is positive where displacements in y directions (v) is 
negative. The value of displacement of Pd matrix is always higher than the AI 2 O 3 due 
higher coefficient of thermal expansion (CTEs) of Pd ( a m = 1 1.5 x 10 " 6 /K) than that of 
AI 2 O 3 (Oc = 8 x 1CT 6 /K). In 01 case, all the displacements values are negative it means 
that the particles are moving in negative x and y directions. Displacements are increasing 
with the increasing of distance. 

Figures 4. 6 . a and 4.6.b show the variation of displacement for nodal points at y- 
axis line for both 10 and 01 cases respectively at AT = 800K. The displacements value 
for initial part does not change much but at the interface the value changes. The Pd 
displacements are higher due to high CTEs of Pd matrix. 


b 



Temperature in Kelvin 



Figure 4.3 — Temperature profile as a function of radial 
distance at AT= 800K. 


S7 



m A1& 



Figure 4.4. a - Cylindrical Precipitate is inside and 
Pd matrix is outside at AT=800K (Ceramic volume 
fraction 0.25) i.e. 10 Case 



Radial Distance (Mts) 


Figure 4.4.b - Cylindrical Precipitate is outside and 
Pd matrix is inside at AT=800K (Ceramic volume 
fraction 0.75) i.e. 01 Case 
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Figure 4.5.a -- Variation of displacements at nodal 

points as a function of distance along x-axis at 
AT=800K for IO case. 
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displacements in metres 



Figure 4.5.b — Variation of displacements at nodal 

points as a function of distance along x-axis at 
AT=800K for OI case. 
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Figure 4.6.a -- Variation of displacements at nodal 

points as a function of distance along y-axis at 
AT=800K for 10 case. 
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Figure 4.6.b — Variation of displacements at nodal 

points as a function of distance along y-axis at 

AT=800K for 01 case. 
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4.4 Estimation of Strains 


After having evaluating displacements in x and y directions, strains can be 
evaluated from equation 4.17. Strain {s} e matrix at each element is evaluated at average 
gauss points. We have used 4 gauss points to integrate the [k] e and {f *} 2 matrix. Average 
values of gauss points are made and these values were used to evaluate all the results 
involve elemental composite. 

Figure 4.7 shows the variation of £xx as a functions of x-axis distance for 10 and 
01 cases. The magnitude of strains for Pd matrix is higher than AI 2 O 3 . The values of 
strains at 2 nd , 5 th & 8 th divisions are nearly same for both the cases. Strains decreases with 
increase of distance because the temperature is also decreasing with increase of distance 
(equation 4.2). the lesser the temperature within the composite lesser will be the strain 
developed. Fluctuation of 2 nd division at interface is more as compared to other division. 

Figure 4.8 shows the variation of £yy as a distance for 10 and 01 cases at AT = 
800K. The trend obtained for Figure 4.7 nearly same as Figure 4.6 but the values are 
different. The strains for Pd is always more then AI 2 O 3 due to high coefficient of thermal 
expansion of AI 2 O 3 . 

Figure 4.9 shows strains (e xy ) as a function of distance at x-axis for both the cases 
at AT = 800K. The values of 01 case are higher than 10 case. CTEs of Pd (a m = 1 1.5 x 
10 ' 6 /K.) is higher than AI2O3 (a c = 8 x 10' 6 /K), Pd’s strains are always high as compared 
to AI 2 O 3 . 


4.5 Estimation of Stresses 

Stresses can easily be evaluated if the strains and temperature profile at each 
element are known. It is obtained from equation 4.8. Stresses depend on strains, 
temperature and material properties per element. All the figures are evaluated at average 
of gauss points. 
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Figure 4.7 -- Variation of strain (8 n ) as a function of 
distance along x-axis for AT=800K for IO and OI 


cases. 



Strain s 



Figure 4.8 — Variation of strain (Syy) as a function of 
distance along x-axis for AT=800K for IO and OI 


cases. 




strain 



Figure 4.9 - Variation of strain (e^) as a function of 
distance along x-axis for AT=800K for IO and OI 
cases. 
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Figure 4.10. a and Figure 4.10.b show the variation ct xx curve as function of 
distance along x-axis for both 10 and 01 cases at AT = 800K. For 10 case, initially as the 
temperature decreases with distance, stresses start decreasing but as the a xx reaches the 
Pd-A^Ch interface, AI2O3 precipitate stresses shoot up due to high Young’s Modulus (E c 
= 380 GPa) as compared to Pd matrix (Em = 123 GPa). To compensate the tensile stress 
in AI2O3, stresses developed in Pd matrix is compressive most of the time, this is shown 
by elements at the 5 th div. For 10 case, the stresses obtained by most of AFO3 are tensile, 
whereas the stresses for Pd matrix is compressive. The reason behind this is that A^CVs 
Young’s modulus value is much higher than Pd matrix. 

Figure 4.11. a and Figure 4.1 1 .b show the variation a xx curve as function of 
distance along x-axis for both 10 and 01 cases at AT = 800K. The values obtained are 
similar to Figure 4. 9. a and 4.9.b. The stresses for element of 2 nd division fluctuates most. 

Figure 4.12 shows the <j xy stress curve as the function of x-axis at AT = 800K. The 
results obtained for 2 nd , 5 th and 8 th divisions are nearly same. The magnitude of a xy is 
higher by a factor of 10 as compared to o H and ct > 7 . Stresses (a xy ) of AI 2 O 3 is high due to 
high Young’s Modulus of A1 2 0 3 (E c = 380 GPa) than that of Pd matrix (E m = 123 GPa). 
The decreasing trend of stresses occurred due to decreasing of temperature with distance. 

4.5.1 Estimation Of Polar Stresses From Cartesian Stresses 

Polar stresses are important for calculating solubility enhancement ratio as it 
involves redial stress. Polar stresses are evaluated from the cartesian stresses from the 
following equations 


= 

cos 2 9 + < 7 VV sin 2 9 + 2a ^ sin<9cosi9 

4.27 

= <3\,v 

sin 2 9 + cos 2 9 -2a xy sin 9 cos 6 

4.28 

= (CT 

n , - cr XK ) sin 9 cos 9 + a n , (cos 29) 

4.29 
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10 case 



Figure 4. lO.a — Variation of stress (a n ) curve as a 
function of distance along x-axis at AT=800K for IO 


case. 




stress (MPa) 



Figure 4.10.b — Variation of stress (a n ) curve as a 
function of distance along x-axis at AT=800K for OI 


case. 




stress (MPa) 



Figure 4.1 l.a - Variation of stress (Gyy) curve as a 
function of distance along x-axis at AT=800K for IO 


case. 





Figure 4.11.b — Variation of stress (ciyy) curve as a 
function of distance along x-axis at AT=800K for OI 
case. 
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stress (M Pa) 



Figure 4.12 - Variation of stress (a^y) curve as a 
function of distance along x-axis at AT=800K for IO and 
01 cases. 
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The radial stress is used for calculating tensile hydrostatic stress for understanding 
solubility enhancement ratio. Radial stress (cin-) means the stress is evaluated at ‘r’ 
distance whereas tangential stress (a 00 ) at the ‘0’ direction. Figure 4.13.a and Figure 
4.13.b shows the variation of radial stress with radial distance for 10 and 01 cases 
respectively at AT = 800K. 

The magnitude of radial stress of AI 2 O 3 is more compared to Pd due to high 
Young’s modulus of AI 2 O 3 . Average value of gaussian values are taken, the results 
obtained by 5 th , 6 th , 7 th and 8 th divisions are same compared to 4 th , 3 rd , 2 nd and 1 st 
divisions. So, results of 1 st , 2 nd , 3 rd and 4 th division are plotted in Figure 4.13. a and 4.13.b. 

In both 10 and 01 case, the major radial stresses are generated at elements at 3 rd 
and 4 th divisions and very little is contributed from 1 st and 2 nd divisions. The value of 
radial stresses for 10 and 01 cases become more negative to less negative as the distance 
increases from origin, this is due to decrease of temperature from inside of composite to 
outside of composite. The stresses obtained for 10 case is more negative than obtained 
for 10 case. 

Figure 4.14.a and Figure 4.14.b show tangential stress (aee) curve as a function of 
radial distance at AT = 800K for 10 and 01 cases. The tangential stresses depend on ‘0’ 
values with respect to x-axis. The figures show symmetry with respect to ‘0’ values. The 
stresses obtained for elements at 1 st , 2 nd , 3 rd and 4 th divisions are symmetry with 8 th , 7 th , 
6 th and 5* divisions respectively. The results obtained for elements at 1 st to 4 th division 
are positive where as the values for elements at 5 to 8 divisions are negative. The 
magnitude of tangential stresses for elements at 1 st and 8 th division are more than stresses 
for 4 th or 5 th division. The results obtained for 10 case is more compared to 01 case, this 
is due to previous values of cr xx , oyy and Oxy- The stresses generated at AI 2 O 3 precipitate is 
more as compared to Pd matrix and it is due to large value of Young’s modulus of AI 2 O 3 . 



4.6 Estimation Of Solubility Ratio Curve 

The tensile hydrostatic stresses are generated from radial stresses as given in 
equation 3.8. Tensile hydrostatic stresses obtained in Pd matrix is averaged and averaged 
hydrostatic stress value is put in the equation 1 . 1 to obtain the solubility ratio. The ratio 
becomes one for 100% AI2O3 or 100% Pd matrix, the radial stresses generated for 100% 
particle is zero. The above ratios are obtained for particular ceramic volume fraction 
keeping temperature difference constant. Ratios are obtained by varying ceramic volume 
fraction and these results are plotted in Figure 4.14.a and Figure 4.l4.b for both IO and 
OI cases respectively. 

Figure 4. 15. a and Fig.4. 15. b shows ratio as 1 for 100% precipitate or 100% 
matrix. The ratio increases with the higher temperature difference. In IO case the ratio 
slowly decreases due to the decrease of Pd percentage in the composite. 

4.7 Estimation of Helmholtz Free Energy 

To calculate the potential in composite, Helmholtz free energy density in plane 
strain condition (eq. 2.7) needs to be written in the component form. 


V = 


E 

2(1 +u) 




1 Eo / 

2 (l+u)(l-2u) Sjx 





4.28 


Helmholtz free energy over the total volume is given by 

\\f = ^.dA 4.29 

D 

As the radial distance is taken for 0 to 12mm, so the thickness is taken as 1mm. 

At particular A1 2 0 3 volume fraction, HFEs are obtained for fixed tempertaure difference 
and these results are plotted in Figures 4. 16. a and 4.16.b. Figure 4.16. a and Figure 4.16.b 
show the variation of HFEs as a function of ceramic volume fraction for IO and OI cases 
at different temperature difference values. 


radial stress (MPa) 


10 case 



Figure 4.13.a -- Variation of radial stress (an-) curve 
as a function of distance along x-axis at AT=800K for 
IO case. 
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Figure 4.13.b -- Variation of radial stress (an-) curve 
as a function of distance along x-axis at AT=800K for 
OI case. 





Figure 4.14.a — Variation of tangential stress (aee) curve 
as a function of distance along x-axis at AT=800K for 
10 case. 



tangential stress (GPa) 



0.000 0.002 0.004 0.006 0.008 0.010 0.012 0.014 

Radial distance along x-axis (Metres) 


Figure 4.14.b - Variation of tangential stress (aee) curve 
as a function of distance along x-axis at AT=800K for 
OI case. 





4.6 Estimation Of Solubility Ratio Curve 

The tensile hydrostatic stresses are generated from radial stresses as given in 
equation 3.8. Tensile hydrostatic stresses obtained in Pd matrix is averaged and averaged 
hydrostatic stress value is put in the equation 1.1 to obtain the solubility ratio. The ratio 
becomes one for 100% AI2O3 or 100% Pd matrix, the radial stresses generated for 100% 
particle is zero. The above ratios are obtained for particular ceramic volume fraction 
keeping temperature difference constant. Ratios are obtained by varying ceramic volume 
fraction and these results are plotted in Figure 4.14.a and Figure 4.14.b for both IO and 
OI cases respectively. 

Figure 4. 15. a and Fig.4.15.b shows ratio as 1 for 100% precipitate or 100% 
matrix. The ratio increases with the higher temperature difference. In IO case the ratio 
slowly decreases due to the decrease of Pd percentage in the composite. 


4.7 Estimation of Helmholtz Free Energy 

To calculate the potential in composite, Helmholtz free energy density in plane 
strain conditions (equation 2.7) needs to be written in the component form 


¥ = 


2(1 + 0 ) 


+ £ 


yy 


' 2 £ xy 2 ) + T 


Ev 


2 (1 + d)(1 - 2v) 


( s xx + £ x>' ) ^ _ 2 V ) + £ >y ) 


4.28 


Helmholtz free energy over the total volume is given by 

V(/ = \y/.dA 4.29 

D 

As the radial distance is taken for 0 to 12mm, so the thickness is taken as 1mm. 

At particular AI2O3 volume fraction, HFEs are obtained for fixed tempertaure difference 
and these results are plotted in Figures 4. 16. a and 4.16.b. Figure 4.16. a and Figure 4.16.b 
show the variation of HFEs as a function of ceramic volume fraction for IO and OT cases 
at different temperature difference values. 


7S 



At 0% A1 2 0 3 precipitate, Pd matrix has 100% HFE. With the increase of A1 2 0 3 
ceramic precipitate, contribution of AI 2 O 3 HFEs slowly increases, where for Pd matrix it 
begins to decrease. At 100% AI 2 O 3 precipitate, no contribution is made by Pd matrix. The 
total HFEs energy is obtained simply by adding energies of Pd matrix and A1 2 0 3 
energies. The total energies is mainly contributed by Pd matrix as compared to AI 2 O 3 
precipitate. 

Energies obtained from IO case is more than OI case due to large values obtained 
from strains for IO cases. Energies depend on the strains and temperature and not on 
stress values. As Pd strains are higher so their HFEs is much higher than A1 2 0 3 
precipitates. 




Figure 4.15.a ~ Variation of solubility enhancement 
curve as a function of ceramic volume fraction for 
different temperature difference values for IO case. 



Solubility Enhancement ratio 



Al 2 0 3 Ceramic volume fraction 


Figure 4.15. b — Variation of solubility enhancement 
curve as a function of ceramic volume fraction for 
different temperature difference values for OI case. 
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Figure 4.16.a — Variation of Helmholtz free energy 
as function of ceramic volume fraction at different 
temperature difference values for IO case. 
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Chapter 5 

CYLINDRICAL PRECIPITATES (CASE II) 

Now we will considered intermediate ceramic precipitates to study the thermal 
residual stresses in composites. Figure 5.1 shows discretization of intermediate 8 
elements precipitates in infinitely long cylindrical matrix. Results are obtained by using 8 
elements AI2O3 precipitates at different locations. Here radical stress, solubility 
enhancement and Helmholtz free energy results are evaluated to study the above figure. 

Figure 5.2 shows the variation of radical stress as a function of radical distance of 
8 elements precipitates. Radical stresses at middle of composite has high magnitude 
value due to presence of AI 2 O 3 precipitates. The stresses slowly increases as we slowly 
proceed from x-axis to a direction 45° to x-axis. Maximum magnitude of stresses is 
shown. 

Different ratios are obtained at various ceramic volume fraction at constant 
temperature difference. These results are plotted in Figure 5.3 for different temperature 
differences. The ratio has value 1 for 0% AI 2 O 3 ceramic fraction. As the ceramic volume 
fraction is increased slowly, the changes in ratio is very small. High temperature 
difference results in high residual stresses which results in high solubility ratio. 

Figure 5.4 shows HFEs per mole as a variation of ceramic volume fraction at 
different temperature difference. The ceramic volume fraction only changes from 0 to 
0.15, contribution of AI 2 O 3 precipitate on HFEs as compared to Pd matrix is nil. So, only 
HFEs of Pd matrix and total composites are plotted. HFEs are calculated using equation 
(2.7). High energies are produced for high temperature difference due to high thermal 
strains in composites. 
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Figure 5.1 - Discretization of the quarter domain using 
2-D elements having 8-elements intermediate 
precipitates. 
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Figure 5.2 -- Variation of radial stress (cr rr ) as a function 
of distance along x-axis at AT=800K for case II. 
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Figure 5.3 -- Variation of solubility enhancement curve 
as a function of ceramic volume fraction for different 
temperature difference values for case II. 
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Figure 5.4 -- Variation of Helmholtz free energy as 
function ceramic volume fraction at different 
temperature difference values for case II. 
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CYLINDRICAL PRECIPITATES (CASE III) 

Only 4 AI 2 O 3 precipitates are embedded at the middle infinitely long composite at 
the middle infinitely long composite. Figure 6 . 1 shows discretization of composite with 
intermediate 4 ceramic precipitates elements. Here, whole matrix is influenced by only 4 
AI 2 O 3 precipitates. 

Figure 6.2 shows the variation of radial stress as the function of radial distance at 
AT = 800K. AI 2 O 3 precipitate contains in 3 rd and 4 th divisions, radial stresses at 3 rd and 4 th 
divisions are more negative. No changes occur for elements at 1 st and 2 nd division as it 
contains Pd matrix only. 

Figure 6.3 shows solubility of intermediate 4 elements precipitates as a function 
of radial distance at different temperature difference. The ratio increases with more 
temperature difference. The ceramic volume fraction is less compared to precipitates 
Figure 5.1. The solubility ratio in present case is much more than case II and is shown in 
Table 6.1. 

The volume fraction of Pd matrix of case IH is more than that of case n, so HFEs 
of case III is higher than case II. The contribution of HFEs is very small because its 
volume fraction is less than 1%. HFEs graph is plotted on Figure 6.4. As the percentage 
of AI 2 O 3 is very less, helmoltz free energy is contributed mainly from Pd matix. In 
Figure 6.4, HFEs of precipitates are neglected due to low value. 



Radial distance 



Figure 6.1 - Discretization of the quarter domain using 
2-D elements having 4-elements intermediate 
precipitates. 
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Figure 6.2 — Variation of radial stress (a rr ) as a function 
of distance along x-axis at AT=800K for case EDE. 
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Solubility Enhancement 



Fig.™ “ - "'solubility enhancement cnrv, 

ns a function of ceramic volume fraction at different 
temperature difference values for case III. 
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Table 6.1 — Comparison of solubility curve for case II 


and case III 
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Figure 6.4 — Variation of Helmholtz free energy as a 
function of ceramic volume fraction at difference 
temperature different values for case III. 





6.1 Limitations of case I, case II and case III. 

It involves plain strain condition so it becomes infinite cylindrical matrix. We 
cannot predict results for 3-dimentional rather it predicts for 2-dimentional. We assume 
that strain at z direction is zero so we basically can’t predict anything about strain at z 
direction. Strain energy evaluated at per unit length and it cannot be evaluated over total 
volume. Elasto-plastic or plastic model are not evaluated to study the plastic behaviour of 
composites. 
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CONCLUSIONS AND SCOPE FOR FUTURE 

WORK 

7.1 Conclusions 

The hydrogen solubility enhancements in Pd-AhCb composites (produced by internal 
oxidation of Pd-Al alloys) due to thermal residual stresses has been evaluated by the 
finite element method. The following are the salient conclusions of the study: 

1. Two different types of systems were evaluated: spherical and cylindrical. The FEM 
approach was applied to both the methods. 

2. In the case of the spherical system, the alumina was assumed to be within the Pd 
matrix (10 case) as well as surrounding the Pd phase (01 case). Results were obtained 
for both the cases. No temperature variation was considered within the Pd-Al 20 3 and 
the temperature difference applied to the difference between high and room 
temperature. 

3. The salient conclusions using the spherical system are: 

(a) For both 10 and 01 cases, displacements for the Pd matrix were always higher 
than those in AI 2 O 3 . This has been explained by the coefficient of the thermal 
expansion of Pd when compare to that of AI 2 O 3 . 

(b) Radial and tangential strains were larger in Pd. This has again been explained by 
the thermal expansion coefficient of Pd. Moreover, radial strains at the Pd-Al 20 3 
interface showed a irregular discontinuity due to change of material properties. 
However, the tangential strains a regular discontinuity . 

(c) The radial stresses developed in the Pd and A1 2 0 3 were compressive for both 10 

and 01 cases. 

(d) Solubility enhancement ratio curve determined for the different volume fraction 
of alumina in the composite. It was found that the ratio exhibited a maximum for 
the 10 case which appeared at relatively low volume fraction of alumina. This has 
been explained by the high constraints in the Pd-Al 2 0 3 phases. Moreover the 





solubility enhancement ratio were higher for the larger AT and this has been 
explained by the higher stresses generated due to large temperature difference. 

(e) Helmholtz free energy (HFE) depends on the strains and the temperature. As the 
strains obtained in the Pd matrix was always higher, this result in higher HFEs for 
Pd matrix. In both 10 and 01 cases, the overall HFEs decreased with increasing 
alumina fraction. 

4. In the case of infinitely long cylindrical system, three types of situations were again 
considered. In the first situation, the alumina precipitates were considered to be 
surrounded by the Pd matrix and the volume fraction of alumina was varied. In the 
second situation, the cylindrical alumina precipitates were placed within the syatem 
such that Pd matrix enclosed it on both the sides. In the third situation, the alumina 
precipitates were placed in specific area of the system. The temperature change 
occurred within the composite. The temperature change within the composite was 
considered by ensuring the temperature in the center and outer surface of the 
composite was different and this was represented by AT. 

5. Salient features for Case I system were as follows: 

(a) The Pd phase displacements were always high due to its high coefficient of 
thermal expansion. Displacements varied continuously in the entire composite due 
to the imposed temperature gradient 

(b) Strain in the Pd were always higher in the 10 case when compared to 01 case 

(c) Stresses (cartesian) generated in the AI 2 O 3 were tensile in nature where as for Pd 
it was compressive. This has been explained based on the material properties and 
matches with previously reported results in the literature. The ctxy values were 
much higher than compared with Oxx and Gyy in both 10 and 01 cases. The 
tangential stresses were symmetric with respect to 0 values. The radial and 
tangential stresses of AI 2 O 3 were always higher than Pd. 

(d) Solubility ratio were higher for larger temperature difference. These were higher 
in 10 case when compared with the 10 case for all volume fraction of alumina. 

(e) The Helmholtz free energy of Pd matrix was higher. Higher temperature 
difference contributed higher HFEs because of enhanced generation if residual 


stresses. 



6. The salient features of the cylindrical system for situation two and three were as 
follows: The ceramic volume fraction was relatively low in these two cases. As the 
solubility enhancement ratios obtained for the situations were higher compared to 
situation one due to large constrained offered by the Pd matrix. Helmholtz free 
energies depend on strains and temperature. The strains for Pd matrix are always 
higher than AI2O3. For Case II and III, AI2O3 volume fraction is very low so the HFEs 
for AI2O3 is very less compared to Pd matrix. Overall HFEs is also higher than Case I 
due to greater Pd volume fraction. 

7. Finite element method is a powerful technique to analyze residual stresses in 
composite systems and the present study enhances this viewpoint. 


7.2 Scope for future work 

1. From the model the strain Szz could not be evaluated. Future work should address this 
aspect. 

2. The models considered here are purely elastic cases. In future work, elasto-plastic and 
plastic models may be evaluated specially in the case of metal matrix composites. 
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